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Abstract 

We estimated the galactic model parameters for a set of 20 intermediate-latitude 
fields with galactic longitudes 0° < I < 100° and 160° < I < 240°, included in the 
currently available Data Release 6 (DR6) of the Sloan Digital Sky Survey (SDSS), 
to explore their possible variation with galactic longitude. The local space densities 
of the thick disc and halo are almost the same for all fields, ((^/rai)} = 6.52% 
and ((713/ni)) = 0.35%, respectively, a result different than the one cited for high- 
latitude fields. The thin disc's scaleheight is 325 pc in the galactic centre changes 
to 369 pc in the third quadrant, which confirms the existence of disc flare, whereas 
the thick disc scaleheight is as large as 952 pc at galactic longitude I = 20° and 10% 
lower at I = 160°, which confirms the existence of the disc long bar in the direction 
I ~ 27°. Finally, the variation of the axis ratio of the halo with galactic longitude 
is almost flat, {(c/a)) = 0.56, except a slight minimum and a small maximum 
in the second and third quadrants, respectively, indicating an effect of the long 
bar which seems plausible for a shallow halo. We estimated the metallicities of 
unevolved G-type stars and discussed the metallicity gradient for different vertical 
distances. The metallicity gradient is d[M/ H]/dz ~ —0.30 dex kpc -1 for short 
distances, confirming the formation of this region of the Galaxy by dissipational 
collapse. However, its change is steeper in the transition regions of different galactic 
components. The metallicity gradient is almost zero for inner halo (5 < z < 10 kpc), 
indicating a formation of merger or accretion of numerous fragments such as dwarf 
galaxies. 

98.35.Ln Stellar content and populations; morphology and overall structure, 98.35.Gi 
Galactic halo 
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1 Introduction 



Our knowledge of the structure of the Galaxy, as deduced from star count data 
with colour information, is more precise nowadays with the advent of new 
surveys. Researchers have used different methods to determine the galactic 
model parameters (cf. iBilir. Karaali fc Gilmord l2006al . Table 1). One can see 
that there is an improvement for numerical values of the model parameters. 
The local space density and the scaleheight of the thick disc can be given as an 
example. The scaleheigh t of the thick disc has s teadily changed o ver the years 
from 1.45 to 0.65 kpc flGilmore fc Reid 19831 ; IChen et all l200ll ) and higher 
local densities (2 - 10%). lJuric et al.l (120081 ) studied SDSS data and obtained 
12% relative local space density and 900 pc scaleheight for the thick disc. The 
authors assumed that 35% of the sample stars were binaries. 



Chen et al. 


( 


2001). 


Siegel et al. 


( 


2002) 



of the thick disc 6.5 - 13% and 6 - 10%, respectively. However, one expects 
the most evolved numerical values from these recent works. That is, either 
the range for this parameter should be small or a single value with a small 
error should be given for it. It seems that researchers have not been able to 
choose the most appropriate procedures for this topic. Large range or differ- 
ent numerical values for a specific galactic model parameter as estimated by 
different researchers may be due to several rea sons: (a) The galactic model 



parameters are absolute magnitude d ependent (IKaraali. Bilir fc Hamzaoglu 



2004 ; IBilir. Karaali fc Gilmord l2006al ). Hence, any procedure which excludes 
this argument galactic model parameters spread in a large range, (b) The 
galactic model parameters are ga lactic latitude /longitude dependent. The two 
works of iBuser et al.l (119981 . Il999l ) confirm this suggestion. Although these au- 
thors give a mean value for each parameter, the values o f a given paramet er 
are not equal for different fields. Several re c ent studies ( Bilir et al. 2006bl lcl; 
Cabrera-Lavers et al.l 120071 ; lAk et al.l l2007al ; IBilir et al.l l2008al ) also support 
this argument, (c) Galactic model parameters change with limiting distance of 
completeness. That is, a specific model parameter is not the sa me for each set 



of ga lactic model parameters estimated for different volumes (IKaraali et al. 
2007t l . 



The difference between disc galactic model parameters estimated for fields with 
different galactic latitudes and longitudes are due to the influence of the disc 
flaring and warping. The disc of the Galaxy is far from being radially smooth 
and uniform. On the contrary, its overall shape presents strong asymmetries. 
While the warp bends the galactic plane upwards in the first and second 
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galactic longitude quadrants (0 < / < 180°) and downwards in the third and 
fourth quadrants (180 < / < 360°), the flare changes the scaleheight as a 
function of radial distance. 



This warp is present in all g alactic components: dust, gas and stars (lLopez-Corredoira et al 



20021 ; iMomany et al.l 120061 ) . The stellar and gaseous flarings for the Milky 



Way are also consistent with each other (IMomany et al.l . 120061 ) , showing that 



the scaleheight of the thin disc (ifi) increases with the galactocentric ra- 



dius for R > 5 kpc dKent. Dame &: Fazio 



Narayan &: Jodl2002l ; lLopez-Corredoira et al 



1991 



2002 



Drimmel & Spereel 


2001: 


Momanv et al. 


2006). The 



behaviour of this flare in the central discs of spiral galaxies has not been stud- 
ied so well due to inherent di fficulties in separating the seve ral contributions 
to the observed counts or flux. lLopez-Corredoira et al.l (120041 ) found that there 
is a deficit of stars compared to the predictions of a pure exponential law in 
the inner 4 kpc of the Milky Way, which could be explained as a flare, which 
displaces the stars to greater heights above the plane as the distance to the 
galactic centre decreases. 

In this scenario, the mean disc (-2=0) can be displaced as much as 2 kpc 
b etwe en the location of the maximum and the minimum amplitudes of the 



warp (IDrimmel fc Spergelll200ll : lLopez-Corredoira et al.ll2002t IMomany et al. 



20061 ). On the other hand, the scaleheight of the stars can sho w differences 
up to 50% of the value for H-\ in the range 5 < R < 10 kpc (lAlvesi l2000l : 
Lopez- Corredoira et al.ll2002l : IMomany et al.ll2006l ) to fit a global galactic disc 
model that accounts for all these inhomogeneities is, at the very least, tricky. 
Because of this, the results in the galactic model parameters might depend on 
the sample of galactic coordinates used, as the combined effect of the warp 
and flare will be different in every direction in the Galaxy and hence along 
different lines of sight. 



There is an additional reason for the difference between the numerical values 
of a given galactic model parameter estimated in different directions of the 
Galaxy, mainly at greater galactocentric distances. These are the observed 
overdense regions with respect to an axisymmetric halo, for which two com- 
peting scenarios h ave been proposed: the fi r st one is concerned with the triax - 
iality of the halo (iNewberg fc Yannv! [20061 : IXu et al.l [20061 : ljuric et alJl2008t l. 
wherea s the second one is relat ed to the remnants of some historical merger 
events ( IWyse &: Gilmord . 120051 ) . 



The size of the field used for galactic model estimation changes from author to 
author. There are star fi elds with sizes ranging from a few square degrees (cf. 
Buser et al.l Il998l . |l999|) to fields with sizes larger than one thousand squ are 
degrees. For example the size of the field investigated by lJuric et al.l (120081 ) is 
6500 deg 2 . Small sizes may involve some gaps causing unreliable galactic model 
parameters. On the other hand, very large fields may be contaminated by rem- 
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nants of mergers. Also, as the galactic model parameters are galactic latitude 
and galactic longitude dependent, a specific model parameter estimated from 
investigation of a large field can be assumed as the mean of corresponding 
model parameters estimated for many neighbour small fields. 

Most of the SDSS data used in the works claimed above for galactic model 
estimation are reduced for high latitude fields. In this paper, we dealt with 
the data of intermediate latitude fields. We fixed the size of the fields with 10 
deg 2 , thus excluding the effect of the field size on the galactic models. Also, 
we avoided the contamination of mergers' remnants. So, we used recent SDSS 
Data Release 6 (DR6) of 20 intermediate latitude fields of equal size with 
galactic longitudes 0° < I < 100° and 160° < I < 240° and estimated galactic 
model parameters for thin and thick discs and halo and compared them with 
the ones estimated for high latitude fields. Such a procedure gave us the chance 
to compare the trends of the galactic models estimated at different galactic 
latitudes and longitudes, as well as their interpretations. 

In Section 2, data and reductions are presented. The galactic model parameters 
and their dependence on galactic longitude and metallicity are given in Section 
3. Section 4 provides a summary and discussion. 



2 Data and Reductions 



The SDSS is a large, international collaboration project set up to survey 10 
000 square-degrees of the sky in five optical passbands and to obtain spectra of 
one million galaxies, 100 000 quasars, and tens of thousands of galactic stars. 
The data are being taken with a dedicated 2.5-m telescope located at Apache 
Point Observatory (APO), New Mexico. The telescope has two instruments: a 
CCD camera with 30 CCD chips, each with a resolution of 2048x2048 pixels, 
totaling approximately 120 million pixels, in the focal plane and two 320 fiber 
double spectrographs. The imaging data are tied to a network of brighter 
astrometric standards through a set of 22 smaller CCDs in the focal plane 
of the imaging camera. A 0.5-m telescope at APO has been used to tie the 
imaging data to brighter photometric standards. 

The SDSS obtains images almost simultaneously in five bro ad bands (u, q 
r, i a nd zj}\ centred at 3551, 4686, 6166, 7480 and 8932 A flFukugita et al. 



19961 ) up to an apparent magnitude of 24 .4, 25.3, 25.1, 24.4 and 22 9, re - 



spectively, with a signal to noise ratio of 5 (lAdelman-McCarthy et all 120081 ) 



The data becomes saturated at about 14 mag in g, r and i and about 12 



1 Magnitudes in this paper are quoted in the ugriz system to differentiate them 
from the former ugriz system. 
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mag in u and z. The imaging data are automatically processed through a 
series of software pipelines which find and measure objects and provide pho- 
tometric and astrometric calibrations to produce a catalogue of objects with 
calibrate d magnitudes, positi ons and structure information. The photometric 
pipeline ( ILupton et al.l . l200ll ) detects the objects, matches the data from the 
five filters, and measures instrumental fluxes, positions, and shape parame- 
ters (which allows the classification of objects as "point source" , "compatible 
with the point spread function" or "extended"). The SDSS DR6 data can be 
obtained by using the SQL interface and includes the complete imaging of 
the Northern Galactic Cap. The catalogue contains images and parameters 
of 287 million objects over 9 583 deg 2 , and 1.27 million spectra of sources 
(star, galaxy, quasar, etc.) over 7 425 deg 2 including scans over a large range 
of galactic latitudes and longitudes. The photometric calibration has been im- 
proved with uncertainti es of 2% in u and 1% in q, r, i sig nificantly better than 
previous data releases (lAdelman-McCarthy et all 120081 ) . 



2.1 The Sample 



The data used in this work were taken from the SDSS DR6 WEB serve for 
22 intermediate-latitude fields (44°. 3 < b < 45°. 7) covering different galactic 
longitude intervals (0° < I < 250°). Intermediate star fields with longitudes 
larger than I = 250° were not available due to the SDSS observing strategy. 
SDSS magnitudes u, g, r, i and z were used for a total of 2xl0 6 stars in 22 
fields. Although the fields are equal in size (10 deg 2 ), their surface densities 
(number of stars per square-degree) are not the same, following a specific 
trend with galactic longitude (Fig. 1). This is the first clue of galactic model 
parameters depending on galactic longitude. Owing to the SDSS observing 
strategy, data of stars brighter than g = 14 mag, dereddened apparent g 
magnitude, are saturated, and star counts are not complete for magnitudes 
fainter than go — 22.2 mag. Hence, our work was restricted to the magnitude 
range 15 < go < 22 for the evaluation of the galactic model parameters. 



2.1.1 De-reddening of the Magnitudes 



The E(B — V) colour excess was evaluated individually fo r each subsam- 
ple source by using the maps of ISchlegel. Finkbeiner fc David (119981 ) through 
SDSS query server, E(B — V) was reduced to total absorption A v : 



A v — 3.1E(B - V). 



http:/ /www. sdss.org/dr6/ access/index. html 
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Fig. 1. Star counts at (b) = +45° for 22 fields available in DR6. 

In order to d etermine to tal absorptions for the SDSS bands (A m ), A m /A v 
data given byS jl999h . i.e. 1.593, 1.199, 0.858, 0.639 and 0.459 for m = u, 
g, r, i and z, respectively, were used. Thus, the de-reddened magnitudes, with 
subscript 0, are 



w 

9o = 9 
tq = r 
i = i 

z = z 



u — A, 



A 

j4y 



a- 



Ai, 

A z . 



(2) 
(3) 
(4) 
(5) 
(6) 



The total absorptions A m are available in the SDSS query server. 



All the colours and magnitudes mentioned hereafter will be de-reddened ones. 
Given that the location of the vast majority of our targets are at distances 
larger than 0.4 kpc, it see ms appropriate to apply the full extinction from the 
maps flBilir et all l2008bh . Actually, when we combine the distance r=0.5 kpc 
(z r^QA kpc distance f rom the galactic plane) with the scaleheight of the dust, 
( Marshall et all 120061 . if=125 pc), we find that the total extinction is reduced 
to almost 6% of the value galactic plane. 



2.1.2 The Star Sample 



According to IChen et al.l (120011 ). the distribution of stars in the go/{g — r)o 
colour-magnitude diagram (CMD) can be classified as follows: Blue stars in 
the magnitude range 15 < go < 18 are dominated by thick disc stars with 
a turn-off at (g — r) ~ 0.33, while galactic halo stars become significant for 



g > 18, with a turn-off at (g — r) ~ 0.2; red stars ((g 
dominated by thin disc stars at all apparent magnitudes. 



r) > 1.3), are 
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Fig. 2. Distribution of 60 960 p oint sources in th e (r — i)o/(g — r)o two-colour 
diagram. Thick solid line denotes J uric et al. ( 2008l Vs stellar locus. Thin solid and 
dashed lines represent the distance d = ±0.15 and d = ±0.30 mag from the stellar 
locus (Eq. 7 in the test), respectively. Also, absolute magnitudes and spectral types 
are given in the figure, confronted to (g — r)o colours. 



However, the CMDs and the two-colour diagrams for all objects indicate 
that the st ellar distributions are contaminated by extragalactic objects as 
claimed by I Chen et al.l (120011 ) . The star/extragalactic object separation is 
based on the "stellarity param eter" as returned from the SExthactoh rou- 
tines ( Bertin fc Arnoutsl . [l996h . This parameter has a value between (high 
extended) and 1 (point source). The separation works very well to classify 
a point source with a value greater than 0.8. Ne edless to say, this separa- 
tion depends strongly on seeing and sky brightness. lJuric et al.l (120081 ) applied 
an extra procedure in order to remove hot dwarfs, low-redshift quasars, and 
white/red dwarf unresolved binaries from their sample. This procedure con- 
sists of rejecting objects at distances larger than 0.3 mag from the stellar locus, 
i.e. the relation between (g — r) and (r — i) colours: 



( g - r ) = 1.39(1 - exp[-4.90(r - i) 3 - 2.45(r - if 



1.68(r-i) (7) 
-0.05]). 
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Fig. 3. 50/(5 — r )o colour-magnitude diagram for stars in the field with galactic 
longitude / = 60°. 



The procedure of lJuric et al 



Bilir et al.l l2008ah . However, 



(120081 ) works well for high latitude fields (cf. 
it needs to be refined a bit for lower galactic 
latitude fields. We used Hipparcos' logarithmic local space density of stars with 
4 < M g < 10, i.e. D* = 7.49 (see section 2.1.5), as a constraint and revealed 
that the distance d = ±0.15 mag is more appropriate than the distance d = 
±0.30 mag for the fields with galactic latitude (b) = 45°. Fig. 2 gives the 
(9 ~ r )o/( r — i)o two-colour diagram, the loci and the borders of the star 
sample in the star field with longitude I = 60° . The qo/(q — r)o CMP for the 
same field in Fig. 3 confirms the arguments of IChen et al.l (120011 ) . 



The histogram of the apparent magnitudes of the star sample in Fig. 4 shows 
that the sample is complete for stars with apparent magnitudes 15 < g Q < 22. 
Hence, the estimation of the galactic model parameters is based on these data. 
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Fig. 4. Apparent magnitude histogram for point sources (white area) and star sample 
(shaded area) for the field centred at I = 60°. 

2.1.3 Dwarf-Gaint Separation 



We used the selection criteria of iHelmi et al.l (120031 ) to separate the metal- 



poor giants. These authors define the location of the metal-poor giants by the 
following criteria: r < 19 mag, 1.1 < (u — g) < 2.0, 0.3 < (g — r) < 0.8, 
0.1 < Pi < 0.6, \s\> m s + 0.05, where P 1 = 0.910(w - g) + 0.415(# - r) - 
1.28, s = 0.249w + 0.794^o - 0.555r + 0.24 and m s =0.002. When we apply 
these criteria to the star sample we obtain ~ 5 069 metal-poor giants which 
corresponds to 0.6% of the whole sample. The projection of the location of 
the dwarfs and giants on the (g — r)o/ (r — i)o two-colour diagram for the field 
with longitude I = 60° is given in Fig. 5. 

As we omitted the bright stars (go < 15 mag) due to saturation, we do not 
expect a sufficient number of metal-rich stars, affecting our statistics. Thus, 
the total number of dwarfs in the sample is 854 960. 



2.1.4 Absolute Magnitudes and Distances 

Absolute magnitudes were determined using two differ ent procedures. For ab- 
solute magnitudes 4 < M g < 8 we used the procedure of lKaraali. Bilir fc Tuncel 



(200 5. hereafter KBTl). whereas for 8 < M g < 10 we adopted the procedure of 



Bilir. Karaali fc Tuncell (12005. hereafter BKTl ). In the procedure of KBT, the 



absolute magnitude offset from the Hyades main sequence, AM 9 , is given as 



a function of both (g — r) colour and 5 0A3 UV-excess, as follows: 

AMf = C 3 <^43 + C 2 6 2 QA3 + Cl 5 .43 + c , (8) 
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Fig. 5. Metal-poor giant stars (black dots) in the (r- 
for stars in the field centred at I = 60°. 



i)o/(g — r)o two-colour diagram 



where £0.43 is the UV-excess of a star relative to a Hyades star of the colour- 
index (g — r)o = 0.43 which corresponds to 5o.e arid is determined using the 
colour transformation equations between UBV and SDSS photometry (KBT). 
The coefficients Cj (i—0, 1, 2, 3) are functions of (g— r) colour and are adopted 
from the work of KBT, where AM^ is defined as the difference in absolute 
magnitudes of a program star and a Hyades star of the same (g — r)o colour: 



AMg 



Ml - M. 



H 



(9) 



The absolute magnitude of a Hyades star can be evaluated from the Hyades 
se quence, normalize d by KBT. Thi s proce dure is the one u sed in the works 
of kk et ali (l2007aT ). iKaraali et all fl2007f ) and iBilir et all fl2008ah . and has 
two main advantages: 1) there is no need to separate the stars into different 
populations, and 2) the absolute magnitude of a star is determined from its 
UV-excess individually which provides more accurate absolute magnitudes 
compared with the procedure Ll in-situ" , where a specific CMD is used for all 
stars of the same population. When one uses the last two equations (Eqs. 8 
and 9) and the following one (Eq. 10), which provides absolute magnitudes 
for the Hyades stars it gets the absolute magnitude M* of a star: 



M, 



H 



-2.0987(# - r) 2 - 0.0008(m - gf + 0.0842(# - r)(u - g) (10) 



+7.7557(# 



O.1556(u-0) + 1.9714. 
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The procedure of KBT was defined for the colour range 0.09 < (g — r)o < 
0.93, which corresponds to absolute magnitudes 4 < M g < 8. Hence, for the 
absolute magnitudes interval 8 < M g < 10 we used the equation of BKT, 
which provides absolute magnitudes for late-type dwarfs: 



M g = 5.791(0 - r) + 1.242(r - z)„ + 1.412. (11) 



Stars with faint absolute magnitudes are very useful, since they provide space 
densities at short distances relative to the Sun, which combine the space den- 
sities for bright stars at large distances, and the local densities of Hipparcos 



(IJahreiss fc Wielenl . 119971 ). Thus, we have a sample of stars with absolute mag- 
nitudes 4 < M g < 10, which enables us to evaluate space density functions in 
the heliocentric distance interval 0.5 < r < 30 kpc, which corresponds to a 
range of distances of 0.4 < z < 20 kpc from the galactic plane. This interval 
is large enough to estimate a set of galactic model parameters and test their 
changes with galactic longitude. The absolute magnitudes in question and the 
corresponding spectral types (from early F type to early M type) for the locus 
points in the (g — r) / (r — z)o two-colour diagram are shown in Fig. 2 for the 
field centred at I = 60° as an example. The local space density in the absolute 
magnitude interval 4 < M g < 10 is flat and it attributes to a mean value of 
logarithmic space density D* = 7.49. 

In a conical magnitude-limited volume, the distance up to which intrinsically 
bright stars are visible is larger than the distance up to which intrinsically faint 
stars are visible. This causes brighter stars to be statistically overrepresented 
and the derived absolute magnitu des to be too faint. This effect, known as 
Malmquist bias ( Malmquist . 1920t ) . was formalized into the general formula: 



2 d\ogA(g) 



M(g)=M -a 2 * "\ (12) 
dg 



where M(g) is the assumed absolute magnitude, Mq is the absolute magnitude 
calculated for any star using KBT calibration, a is the dispersion of the KBT or 
BKT calibration, and A(g) is the differential counts evaluated at the apparent 
magnitude go of any star. The dispersion in absolute magnitude calibration 
of KBT and BKT is around 0.25 mag, corresponding to an error of 10% 
in photometric distance. We divided the sample into the absolute magnitude 
intervals (4,5], (5,6], (6,7], (7,8], (8,9] and (9,10], and we applied the Malmquist 
bias to stars in each interval separately. This approach provides uniform space 
densities which is essentially the Malmquist bias. Thus, the corrections applied 
to the absolute magnitudes are 0.005, 0.003, 0.007, 0.008, 0.012 and 0.012 for 
the absolute magnitude intervals cited above. The correction of the Malmquist 
bias was applied to SDSS photometric data used in this work. 
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Combination of the absolute magnitude M g and the apparent magnitude go 
of a star gives its distance r relative to the Sun, i.e. 



[(?-M s ]o = 51ogr-5. (13) 



Gilmore. Wyse fc Jonesl (119951 ) quote an error of ~ 0.2 dex in the derivation 
of [M/H] from the UBV photometry for F/G stars which leads to a random 
uncertainty of 20% in the distance estimation. One expects larger distance 
errors for late spectral type stars. However, (u — g) and (g — r) colours are 
more accurate than the (U — B) and (B — V) colours which mitigate this 
excess error for K stars. The distance of the star to the galactic plane can be 
evaluated using its distance r and galactic latitude b: 

z = r sin b. (14) 



2.1.5 Density Functions 

Logarithmic space densities D* = log D + 10 have been evaluated for the com- 
bination of three population components (thin and thick discs and halo), for 
each field where D = N/AV 1>2 ; AF 1>2 = (7r/180) 2 (D/3)(rf - rf); □ denotes 
the size of the field (10 deg 2 ); r\ and r 2 are the lower and upper limiting 
distances of the volume AV 1)2 ; N is the number of stars per unit absolute 
magnitude; r* = \(rf + rD/2] 1 / 3 is the centroid distance of the volume AVi^; 
and z* = r* sin b, b being the galactic latitude of the field centre. The limiting 
distances of completeness rj and zi can be calculated by substituting gi, rj and 
z\ for g, r and z in Eqs. 13 and 14, where g\ is the limiting apparent magnitude 
(15 and 22 mag for bright and faint stars, respectively). The logarithmic den- 
sity function evaluated for the field centred at the galactic longitude / = 60° 
is given in Fig. 6. 



2.1.6 Density Laws 



In th is work we adopted the density laws of Basle group (IBuser et all 11998 



19991 ) . Disc structures are usually parametrized in cylindrical coordinates using 



radial and vertical exponentials: 

Di(x,z)=ni exp(-\z\/Hi) exp(-(x - Ro)/hi), (15) 



where z = z® + r sin b, r is the distance to the object from the Sun, b is the 
galactic latitude, z ^ is the vertical distance of the Sun from the galactic plane 



(IJuric et all 120081 . 24 pc), x is the projection of the galactocentric distance 



on the galactic plane, Rq is the solar distance from the galactic centre (8 kpc; 
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Reidlll993l ). Hi and hi are the scaleheight and scalelength, respectively, and 
n, is the normalized density at the solar radius. The suffix % is 1 for the thin 
disc, whereas 2 for the thick disc. 

The density law for the sphero id component i s para meterized in different 



forms. The most common is the Ide Vaucouleurd (119481 ) spheroid used to de- 



scribe the surface brightness profile of ellipt i cal ga laxies. This law has been 



deprojected into three dimensions by lYoungl (119761 ) as 



D 8 (R) = n s (exp[10.093(l - R/ R ) 1/4 }/ (R/ R ) 7/8 ) (16) 

(1 - [0.08669/(^/^0]), 

where R is the (uncorrected) galactocentric distance in spherical coordinates, 
and n s is the normalized local density. R has to be corrected for the axial ratio 

(c/a), 

R=[x 2 + (z/(c/a)) 2 } 1 / 2 , (17) 
where 



z = r sin b 



x 



[Rl + r 2 cos 2 b - 2R r cos b cos /] 1/2 , (19) 



with r being the distance along the line of sight and, /, b being the galactic 
longitude and latitude respectively, for the field under investigation. 



3 Galactic Model Parameters and Metallicity Distribution 

3.1 Estimation of the Galactic Model Parameters 



We estimated all the galactic model parameters simultaneously, by fitting 
the space density functions derived from the observations (combined for the 
three population components) to a corresponding combination of the adopted 
population-specific analytical density laws. The faintest stars in this work 
provide space densities at short distances from the galactic plane, z ~ 0.4 
kpc. Hence, it is possible to have reliable extrapolation between them and 
the space density of Hipparcos in the solar neighbourhood, D* = 7.49 in 
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Fig. 6. Observed (dots) and evaluated (thick solid line) space density functions 
combined for stars of all three population components. Also, independent space 
density functions of thin disc (thin dashed line), thick disc (solid line) and halo 
(dotted line) for the star field with I = 60° are shown. 



logarithmic form (IJahreiss fc Wielenl . 1 19971 ). corresponding to the mean of the 
local space densities for stars with 4 < M g < 10. 

We used the classical x 2 statistic to estimate galacti c model parameters, which 



is the most commonly method in r ecent studies (Phleps et al.l 2000 



Chen et al 



Bilir et al. 



2001; ISiegel et al.l 120021 : iDu et all 12003 



20061 ; lJuric et al 



2005 



2008 



2008aj). The comparison of the logarithmic density functions de- 



rived from the observations and the analytical density laws is given in Fig. 
6 for the field with galactic longitude I = 60°, as an example. Xmin shows a 
symmetrical distribution, as it can be seen from Fig. 7, which is given as an ex- 
ample. Hence, the errors of the galactic model parameters could be estimated 
by changi ng a given model p arameter until an increase or decrease by 1 was 
achieved (IPhleps et al.l . |2000| ) . Table 1 lists the galactic model parameters for 
20 SDSS intermediate-latitude (44°. 3 < b < 45°. 7) fields, resulting from the 
fits of the analytical density profiles (fields with galactic longitudes I = 150° 
and I = 250° are omitted from the table due to their small sizes, ~4 deg 2 ). 
The columns indicate: galactic longitude (I), scaleheight of thin (Hi) and thick 
discs (H 2 ), local space densities of the thick disc (ri2/ni) and the halo (n 3 /ni) 
relative to the local space density of the thin disc, scalelength of the thin (hi) 
and thick discs (^2); axial ratio of the halo (c/a), reduced chi-square mini- 
mum (Xmin), standard deviation (s) and the corresponding probability. The 
Xmin va l ues are low, whereas the probabilities are rather high confirming the 
reality of the galactic model parameters. 



We also used a different procedure to test any possible degeneracy in the es- 
timation of the model parameters. First, we estimated the local space density 
and the axial ratio for the halo by comparing the logarithmic space density 
function for z > 5 kpc with the analytical density law of the halo (Eq. 16). 
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Table 1 



Estimated galactic model parameters for 20 star fields. 



1 


Hi 


hi 


H 2 


h 2 


(n 2 /ni) 


(c/a) 


(«3/™l) 




s 


Prob 


(dog) 


(PC) 


(kpc) 


(PC) 


(kpc) 


(%) 




(%) 











325 ± 6 


1.00 ± 0.39 


946 ± 26 


4.48 ± 0.39 


6.56 ± 0.63 


0.545 ± 0.01 


0.33 ± 0.02 


7.95 


0.022 


0.997 


10 


332 ± 6 


1.00 ± 0.22 


952 ± 45 


4.21 ± 0.32 


6.55 ± 0.56 


0.560 ± 0.01 


0.33 ± 0.01 


4.36 


0.017 


0.999 


20 


343 ± 7 


1.00 ± 0.43 


946 ± 25 


5.49 ± 0.74 


6.50 ± 0.96 


0.565 ± 0.01 


0.33 ± 0.01 


5.83 


0.019 


0.999 


30 


348 ± 5 


1.00 ± 0.44 


951 ± 20 


4.00 ± 0.44 


6.46 ± 0.39 


0.570 ± 0.01 


0.33 ± 0.01 


7.92 


0.023 


0.997 


40 


352 ± 6 


1.00 ± 0.53 


941 ± 24 


3.67 ± 0.49 


6.49 ± 0.44 


0.569 ± 0.01 


0.33 ± 0.01 


7.23 


0.022 


0.998 


50 


354 ± 9 


1.00 ± 0.69 


934 ± 22 


3.35 ± 0.52 


6.52 ± 0.51 


0.569 ± 0.01 


0.34 ± 0.01 


8.26 


0.023 


0.996 


60 


352 ± 6 


1.05 ± 0.12 


931 ± 22 


3.52 ± 0.94 


6.52 ± 0.49 


0.564 ± 0.01 


0.39 ± 0.01 


7.94 


0.023 


0.997 


70 


353 ± 7 


1.16 ± 0.31 


921 ± 18 


3.11 ± 1.32 


6.56 ± 0.50 


0.551 ± 0.01 


0.36 ± 0.01 


9.31 


0.025 


0.991 


80 


350 ± 7 


1.90 ± 0.77 


902 ± 18 


2.90 ± 1.53 


6.52 ± 0.43 


0.537 ± 0.01 


0.36 ± 0.01 


7.23 


0.021 


0.998 


90 


331 ± 10 


1.01 ± 0.76 


882 ± 24 


2.65 ± 0.78 


6.40 ± 0.45 


0.520 ± 0.01 


0.36 ± 0.01 


7.27 


0.022 


0.998 


100 


330 ± 9 


1.36 ± 0.47 


883 ± 17 


3.90 ± 0.91 


6.47 ± 0.38 


0.510 ± 0.01 


0.36 ± 0.01 


8.01 


0.022 


0.997 


160 


357 ± 8 


1.38 ± 0.12 


860 ±17 


4.15 ± 0.35 


6.46 ± 0.39 


0.489 ± 0.01 


0.33 ± 0.01 


8.10 


0.023 


0.996 


170 


355 ± 9 


1.48 ± 0.14 


863 ± 17 


4.15 ± 0.39 


6.55 ± 0.39 


0.510 ± 0.01 


0.34 ± 0.01 


8.40 


0.023 


0.995 


180 


362 ± 8 


1.66 ± 0.16 


875 ± 32 


4.22 ± 0.56 


6.62 ± 0.54 


0.537 ± 0.01 


0.34 ± 0.02 


7.00 


0.020 


0.999 


190 


363 ± 8 


1.65 ± 0.25 


877 ± 17 


4.50 ± 0.44 


6.62 ± 0.43 


0.551 ± 0.01 


0.35 ± 0.02 


8.29 


0.023 


0.996 


200 


369 ± 9 


1.68 ± 0.17 


904 ± 23 


4.45 ± 0.48 


6.56 ± 0.40 


0.610 ± 0.01 


0.36 ± 0.02 


6.97 


0.021 


0.999 


210 


368 ± 9 


1.68 ± 0.18 


905 ± 20 


4.47 ± 0.50 


6.59 ± 0.40 


0.654 ± 0.02 


0.39 ± 0.02 


8.75 


0.020 


0.994 


220 


365 ± 8 


1.53 ± 0.17 


895 ± 22 


3.80 ± 0.41 


6.55 ± 0.41 


0.651 ± 0.01 


0.37 ± 0.01 


7.74 


0.022 


0.997 


230 


360 ± 9 


1.23 ± 0.13 


899 ± 20 


4.22 ± 0.52 


6.47 ± 0.37 


0.606 ± 0.01 


0.36 ± 0.01 


8.00 


0.023 


0.997 


240 


354 ± 8 


1.09 ± 0.14 


918 ± 19 


4.38 ± 0.60 


6.50 ± 0.39 


0.577 ± 0.01 


0.36 ± 0.02 


8.19 


0.023 


0.996 



Then, we omitted the space density of the halo, estimated by the correspond- 
ing density law, and compared the new density function for z < 5 kpc with 
the combined density laws of thin and thick discs (Fig. 8b). This procedure 
provides galactic model parameters for thin and thick discs independent of 
the model parameters of the halo. The result of the application of two differ- 
ent procedures shows that the corresponding galactic model parameters for a 
specific population are either identical or differ only by a negligible amount. 
Hence, we may argue that there is no degeneracy in the estimated parameters. 

A similar procedure is applied to the thin and thick discs. We estimated the 
local space densities and the scaleheights of thin and thick discs simultaneously 
by using the space density function for 0.3 < z < 5 kpc and compared the 
model parameters of the thick disc with the corresponding ones, which was 
estimated by using space density function for 1.5 < z < 5 kpc, where the thick 
disc is dominant. Since no significant differences could be observed between 
the compared parameters, we concluded that no degeneracy exists between 
the two discs either. 

3.2 The Trend of the Local Space Density of the Thick Disc and Halo 

The local space densities of the thick disc (n 2 /ni) and halo (n 3 /ni), relative 
to the local space density of the thin disc, were plotted versus the galactic 
longitude in Fig. 9. The trend is flat in both panels with exception of two small 
offsets at I = 60° and I = 210° in panel (b). This invariance against galactic 
longitude is not mentioned in the literature, therefore it is new. However, it is 
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Fig. 7. x 2 distribution of estimated Hi (a), H 2 (b) and (c/a) (c) model parameters 
for the star field with / = 60°. 

different than the results obtained for high latitude fields. For example, in our 
recent works ( Karaali et al. 2007I ; Bilir et al. 2008a ) the local space densities 
of thick disc (r^/wi) and halo (713/ni) vary with galactic longitude (see Section 
4). 



3.3 Dependence of the Scaleheights and Scalelengths of the Discs and Axial 
Ratio of the Halo on the Galactic Longitude 



The variation of the scaleheights of the thin and thick discs, Hi and H 2 , 
respectively, with the galactic longitude (I) were plotted in Fig. 10. Hi assumes 
its minimum (325 pc) in the galactic centre direction and it increases with 
longitude in the first quadrant. In the second quadrant, it gives the indication 
of decreasing, however it is not complete. Finally, the third quadrant displays 
a slight variation of Hi with I, where Hi assumes its maximum (369 pc). 



16 



8 




5 10 15 20 




^ I I I I I I I I I I I I 

1 2 3 4 5 6 

z (kpc) 



Fig. 8. Investigation for degeneracy in model parameters of simultaneously estimated 
populations for the star field with I = 60°. The dots and thick line represent observed 
and theoretical space densities, respectively, (a) Using z > 5 kpc data, only the local 
space density and the axial ratio of the halo are estimated independent of discs, (b) 
Data for z < 5 kpc without halo densities are used to estimate the model parameters 
of thin and thick discs. 



The variance of thick disc's scaleheight (H 2 ) with galactic longitude differs 
from variance of thin disc's scaleheight (Hi). H 2 assumes its maximum (952 
pc) at the galactic longitude / = 20° and it decreases within the first quadrant. 
It assumes the minimum value 860 pc at I — 160° and it increases gradually 
in the third quadrant. One can see a similarity in the trends of scaleheight of 
the thick disc (H 2 ) and surface densities of the fields investigated (Fig. 1). 

The numerical values as well as the trends of the scalelengths of two discs are 
also different (Fig. 11). The scalelength of the thin disc (hi) is almost constant 
1 kpc) in the longitude interval 0° < I < 60°, then it increases with high 
errors up to I — 100°. The trend of Hi in the third quadrant resembles a 
parabolic function with a maximum of 1.69 kpc at I = 200°. 

The scalelength of the thick disc (h 2 ) decreases with galactic longitude in the 
first quadrant and it assumes its maximum and minimum, 5.49 and 2.65 kpc 
at / = 20° and I = 90°, respectively. However, it is rather flat in the third 
quadrant (h 2 ~ 4.26 kpc). 
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Fig. 9. Variation of normalized local space densities of thick disc (a) and halo (b) 
as a function of galactic longitude. 
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Fig. 10. Variation of the scaleheights of the thin (a) and thick discs (b) as a function 
of galactic longitude. 

The variation of the axial ratio of the halo (c/a) is given in Fig. 12. The 
trend of (c/a) gives the indication of a smooth variation (0.55 < (c/a) < 0.60) 
except, perhaps, a slight minimum in the second quadrant and a small devi- 
ation towards large values of (c/a), between galactic longitudes / = 200° and 
I = 230°. That is, the halo component of the Galaxy has a disclike structure 
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Fig. 11. Variation of scalelengths of the thin (a) and thick (b) discs as a function of 
galactic longitude. 
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Fig. 12. Variation of axial ratio of the halo as a function of galactic longitude. 

at the intermediate galactic latitudes which is what we were expecting (see 
Section 4). 



3.4 Effect of the Unsolved Binaries 



A high fraction of stars are actually binary systems and being a binary system 
makes stars appear brighter and redder than normally they are. Different frac- 
tional values (defined as /) can be found in the literature. For example, using 
the data in the Gliese catalogue of nearby stars, iBroschd (119641 ) found a value 
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of / = 0.4 due to his simple model for the resolution c riterion. On considering 
the local (within a distance of 10 pc) binary fraction, iReid I (119911 ) concluded 
that the proportion of binaries among "stars" is consistent with a value rang- 
ing from 30 to 50%. When all sy stems in question are binary stars, i.e. / = 1, 
Kroupa. Gilmore &: Toutl (ll99ll ) found that a single mass function provides 
the best representation of a single luminosity fu nction. However, a smaller 
value can not be discarded with high confidence. lHalbwachs I (119861 ) used all 
available data on binary systems and concluded that the proportion of single 
stars among all stellar systems is at most 23% when spectroscopic binaries 
are taken into account. An extensive long term radial velocity study of the 
Hyades cluster reveals that at least 30% of the cluster stars are spectroscopic 
binaries and that essentially all st ars brighter than th e Hyades mainsequence 
stars are actually binary systems (jGriffm et aL| , [1988J ). 



The effect of binary stars were discussed in ( jKroupa. Tout fc Gilmore 111990 



Kroupa. Tout &: Gilmord Il993l . hereafter KTG90 and KTG93, respectively) 
extensively as well as other effects such as metallicity, age, distance etc. KTG93 
adopt the binary fraction / ~ 0.6 — 0.7 as a reasonable value. They give 
the mentioned combined effects as "cosmic scatter" as a function of (V — I) 
colour in the range 0.5 < (V — I) < 4.5. These authors estimate the scatter 
belonging to binaries alone as a = 0.27 mag, if a fraction / = 0.8 of all stars 
are unresolved binary systems. 



The effect of binary stars on the Galactic model parameters were discusse d 
by many authors (cf. ISiegel et al.l 120021 ; lJuric et al.ll2008l ; llvezic et al.ll2008l ). 
The net effect of binarism, stated in the literature, is an underestimation 
of the scaleheights for thin and thick discs. We adopted three fractions of 
binary stars, i.e. 25%, 50% and 75%, and evaluated all the Galactic model 
parameters for thin and thick disc s, and halo f or one of our fields (I = 60°). 
The procedure is similar to that of ISiegel et al.l (120021 ) . For a fraction of stars, 
we added a companion chosen at random to be equal in mass to the primary. 
The colour and magnitude of each primary star was then reset to the combined 
characteristics of the binary. Then, we evaluated the density function of the 
stars using the original calibrations and inferred density law was compared 
with the input value (Fig. 13). The result is a steepening of the measured 
density law and consequent underestimation of the scaleheight. Whereas effect 
of binarism causes an overestimation for the local space densities (Table 2). 



According to our results, a fraction of 50% unresolved binaries decreases the 
scaleheight of the thin disc by 22%, for a spec ific fi eld (I = 60° ) , w hich is 



between the values stated by lJuric et al.l (120081 ) and ISiegel et al.l (120021 ). i.e. 
20 and 25%, respectively. However, the decrease of the scalehe ight for the 
thick disc due to binarism is higher, 45%, than the one stated by ISiegel et al. 
(l2002h . 20-29%. 
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Fig. 13. The effect of unresolved binaries upon the derived density law for different 
ratio of equal mass binaries. The panels show the actual density law (solid line) 
against the derived density law (dashed line) for various binary fraction. 

Table 2 

The effect of unresolved binaries upon the Galactic model parameters for thin and 
thick discs, and halo, for three fractions of binary stars for the field with I = 60° 
(distances in pc). 







f= 


=0.25 




f=0.50 




f= 


=0.75 




Parameter 


Actual Value 


Dcricvcd Va] 


ue Variation 


(%) Dcricvcd Value Var 


ation 


(%) Dcricvcd Val 


uc Variation (%) 


-Hi (pc) 


352 


376 




7 


428 


22 


464 




32 


H 2 (pc) 


931 


1069 




IS 


1350 


45 


1465 




57 


(c/a) 


0.56 


0.55 




2 


0.53 


5 


0.51 




9 


n 3 /m (%) 


6.52 


5.51 




IS 


4.3 


34 


3.9 




40 


n 3 /m (%) 


0.39 


0.38 




3 


0.37 


5 


0.36 




8 



3.5 Metallicity Distribution 



The photometric metallicities of about 130 000 G-type stars with absolute 
magnitudes 5 < M g < 6 were evaluated using the following equation of KBT: 

[M/H] = 0.10 - 3.545 .43 - 39.63(^.43 + 63.51^ 43 . (20) 
This equation was calibrated for the main-sequence stars with 0.12 < (g— r)o < 
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Fig. 14. Metallicity distributions as a function of vertical distance z* for the field 
with I = 60°. 



0.95 which covers the (g — r) colour indices of our sample. Here, £0.43 is the 
normalized UV-excess in SDSS photometry corresponding to 5 . e in the UBV 
photometry. KBT give the range of the metallicity as —2.76 < [M/H] < 0.20 
dex, corresponding to < 5 0A3 < 0.33. 

The metallicity distribution for each field is determined using a mean metal- 
abundance value determined for each of the following distance intervals (in 
kpc): (0.5, 1.5], (1.5, 2], (2, 2.5], (2.5, 3], (3, 4], (4, 5], (5, 6], (6, 7], (7, 8], 
(8, 9], (9, 12], (12, 15]. Fig. 14 shows the metallicity distribution for the field 
with galactic longitude / = 60°, as an example. The projection of the centroid 
distances (r*) of these intervals onto the vertical direction, i.e. z* = r*sin&, 
are as follows: 0.85, 1.26, 1.61, 1.96, 2.52, 3.22, 3.92, 4.62, 5.33, 6.03, 7.57, 
9.66 kpc. The mean metal abundance is adopted as the mode of the Gaussian 
curve fitted to each distribution, except three intervals, (4,5], (5,6] and (6,7], 
where the metallicity distribution is almost flat, and where the median was 
assumed as the mean of the metallicity distribution. The variation of the mean 
metallicities with z* distances for 20 fields is given in Table 3. The metallicity 
histogram, obtained from Table 3, for all fields is given in Fig. 15 for each z* 
distance. 



We coupled the mode/median and the corresponding z* distance for each 
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histogram (Table 4) and plotted them in Fig. 16 in order to investigate the 
vertical metallicity gradient for 20 intermediate latitude fields. One can see 
four different trends in the figure: (1) for relatively short z distances, i.e. 
z* < 2.5 kpc, the variation of [M/H] is rather smooth; (2) for intermedi- 
ate z* distances, i.e. 2.5 < z* < 4 kpc, the variation is steeper but still 
smooth; (3) for 3.5 < z* < 5.5 kpc, the variation is the same as for the 
short distances; (4) and finally for 6 < z* < 10 kpc there is no variation. 
The metallicity gradient describing the first trend, d[M/H]/dz = —0.32 dex 
kpc -1 , is in agreement with the canonical metallicity gradients for the same 
z* distances, and is likely the signature of this galactic regions formation by a 
dissipative collapse. The description of the second trend in terms of a metal- 
licity gradient, d[M/H]/dz = —0.56 dex kpc" 1 , corresponds to the average 
metallicity differences between the three population components involved. The 
third trend, d[M/H]/dz = —0.20 dex kpc -1 , corresponds to the metallicity 
difference between thick disc and inner halo, and finally the fourth trend, 
d[M/H]/dz = —0.01 dex kpc -1 is the very low metallicity gradient of the in- 
ner spheroid. The mean metallicity gradient for the whole z* interval, z* < 10 
kpc, d[M / H\/ dz = —0.25 dex kpc -1 , is i n agreem e nt with the ones appearing 
in the literature for high latitude fields ( Ak et al. . 2007al lbl). 
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Fig. 16. Metallicity versus vertical distance z* . The dashed line represents the 
mean metallicity gradient. 

Table 3 

The metallicity distribution as a function of distance to the galactic plane, for 20 
intermediate latitude ((b) = 45°) fields at different galactic longitudes (0° < / < 100° 
and 160° < / < 240°). Uncertainties (given in parentheses) refer to the last quoted 
digits. Distances are in kpc. 



l/z* 
dcg/kpc 


0.85 


1.26 


1.61 


1.96 


2.52 


3.32 


3.92 


4.62 


5.33 


6.03 


7.57 


9.66 


N 





-0.58(06) 


-0.82(05) 


-0.90(05) 


-0.97(07) 


-1.35(02) 


-1.55(02) 


-1.62(02) 


-1.61(02) 


-2.34(09) 


-2.19(19) 


-2.27(16) 


-2.14(23) 


10058 


10 


-0.34(03) 


-0.58(05) 


-0.68(04) 


-0.70(07) 


-0.83(14) 


-1.50(02) 


-1.64(03) 


-1.80(03) 


-1.95(19) 


-2.20(09) 


-1.98(15) 


-2.34(05) 


13623 


20 


-0.28(06) 


-0.49(09) 


-0.52(05) 


-0.62(05) 


-0.82(19) 


-1.30(02) 


-1.60(03) 


-1.76(03) 


-1.85(09) 


-1.78(13) 


-2.24(09) 


-2.24(11) 


13607 


30 


-0.41(05) 


-0.59(04) 


-0.65(07) 


-0.69(05) 


-0.80(02) 


-1.37(02) 


-1.69(03) 


-1.68(03) 


-2.17(08) 


-2.22(29) 


-2.04(09) 


-2.40(02) 


12765 


40 


-0.54(05) 


-0.65(04) 


-0.75(04) 


-0.78(06) 


-0.81(07) 


-1.44(03) 


-1.70(03) 


-1.76(03) 


-1.74(04) 


-2.09(17) 


-2.15(08) 


-2.36(10) 


11771 


50 


-0.46(05) 


-0.61(05) 


-0.70(05) 


-0.66(05) 


-0.84(09) 


-1.32(03) 


-1.57(03) 


-1.84(04) 


-2.28(09) 


-2.30(10) 


-2.41(43) 


-2.41(09) 


10368 


60 


-0.42(03) 


-0.63(06) 


-0.71(03) 


-0.82(05) 


-0.80(11) 


-1.24(03) 


-1.63(04) 


-1.89(04) 


-2.00(20) 


-2.37(12) 


-2.23(07) 


-2.17(07) 


8430 


70 


-0.47(05) 


-0.52(06) 


-0.61(03) 


-0.75(05) 


-1.00(03) 


-1.37(04) 


-1.64(04) 


-1.74(05) 


-1.77(05) 


-1.76(06) 


-2.27(08) 


-2.10(03) 


7101 


80 


-0.38(04) 


-0.56(05) 


-0.63(04) 


-0.77(04) 


-0.99(03) 


-1.33(04) 


-1.71(05) 


-1.77(05) 


-2.16(06) 


-2.17(25) 


-2.40(16) 


-2.30(14) 


5648 


90 


-0.34(06) 


-0.29(05) 


-0.44(06) 


-0.49(16) 


-0.90(03) 


-1.23(04) 


-1.50(05) 


-1.70(06) 


-1.86(07) 


-2.24(25) 


-2.36(24) 


-2.27(13) 


4102 


100 


-0.32(08) 


-0.51(05) 


-0.59(07) 


-0.74(07) 


-0.79(08) 


-1.30(05) 


-1.66(06) 


-1.89(06) 


-2.00(07) 


-2.21(13) 


-2.10(17) 


-2.30(04) 


3635 


160 


-0.34(07) 


-0.36(12) 


-0.66(12) 


-0.71(12) 


-1.00(04) 


-1.30(06) 


-1.73(07) 


-1.70(07) 


-1.80(07) 


-1.84(19) 


-2.13(19) 


-2.20(17) 


2701 


170 


-0.32(05) 


-0.40(08) 


-0.49(07) 


-0.46(16) 


-0.94(05) 


-1.37(06) 


-1.65(07) 


-1.88(07) 


-2.00(06) 


-2.18(05) 


-2.20(04) 


-2.19(10) 


2492 


180 


-0.38(06) 


-0.43(10) 


-0.42(11) 


-0.64(09) 


-0.99(05) 


-1.33(06) 


-1.70(07) 


-1.80(05) 


-1.86(08) 


-2.19(17) 


-2.16(27) 


-2.30(06) 


2704 


190 


-0.26(09) 


-0.38(05) 


-0.44(20) 


-0.44(09) 


-0.88(05) 


-1.30(06) 


-1.73(06) 


-1.88(07) 


-2.23(09) 


-2.20(33) 


-2.27(28) 


-2.24(21) 


2761 


200 


-0.18(03) 


-0.32(11) 


-0.39(03) 


-0.39(23) 


-0.86(05) 


-1.20(05) 


-1.74(06) 


-1.80(06) 


-1.83(07) 


-1.90(12) 


-2.18(09) 


-2.43(11) 


3062 


210 


-0.30(07) 


-0.48(07) 


-0.60(07) 


-0.62(10) 


-1.05(04) 


-1.31(05) 


-1.79(06) 


-1.96(06) 


-2.20(04) 


-2.24(15) 


-2.29(14) 


-2.41(23) 


3684 


220 


-0.37(07) 


-0.35(09) 


-0.58(06) 


-0.61(08) 


-0.98(04) 


-1.17(06) 


-1.61(06) 


-1.87(07) 


-1.98(06) 


-2.40(07) 


-2.20(03) 


-2.20(04) 


3570 


230 


-0.24(08) 


-0.54(03) 


-0.64(05) 


-0.77(03) 


-0.97(04) 


-1.30(05) 


-1.67(06) 


-1.85(06) 


-2.31(31) 


-2.31(12) 


-2.49(18) 


-2.69(06) 


3526 


240 


-0.33(06) 


-0.52(06) 


-0.58(06) 


-0.61(08) 


-1.04(04) 


-1.35(05) 


-1.70(05) 


-1.93(06) 


-2.00(05) 


-2.22(17) 


-2.28(08) 


-2.30(03) 


3761 



4 Summary and Discussion 



We estimated the galactic model parameters for 20 intermediate-latitude fields 
with galactic longitudes 0° < I < 100° and 160° < I < 240° to explore their 
possible variation with galactic longitude. We evaluated the metal abundance 
for about 130 000 G-type main-sequence stars with absolute magnitudes 5 < 



24 



Table 4 

Mean metallicities (for different distance intervals) calculated from the metallicity 
distributions of 20 fields. Symbols: r: distance from the Sun; r*: centroid distance 
corresponding to the interval r\ — r<i\ z: projection of r* onto the vertical direction; 
[M/H]: mean metallicity. 



ri — T2 


r 


* 

z 


[M/H] 


(kpc) 


(kpc) 


(kpc) 


(dex) 


0.5-1.5 


1.21 


0.85 


-0.35(06) 


1.5-2.0 


1.79 


1.26 


-0.50(07) 


2.0-2.5 


2.28 


1.61 


-0.59(07) 


2.5-3.0 


2.77 


1.96 


-0.69(09) 


3.0-4.0 


3.57 


2.52 


-0.89(06) 


4.0-5.0 


4.55 


3.32 


-1.31(04) 


5.0-6.0 


5.55 


3.92 


-1.67(05) 


6.0-7.0 


6.54 


4.62 


-1.79(05) 


7.0-8.0 


7.53 


5.33 


-1.95(09) 


8.0-9.0 


8.53 


6.03 


-2.22(16) 


9.0-12.0 


10.71 


7.57 


-2.21(15) 


12.0-15.0 


13.66 


9.66 


-2.24(10) 



M g < 6 to investigate the metallicity gradient of different components of the 
Galaxy, as well. 



4-1 Local Space Densities 



The local space densities of the thick disc and halo relative to the local space 
densities of the thin disc with galactic longitude do not vary. This is an in- 
novation and it confirms its definition. However its trend is different than 



the re sults obtained for high latitude fields (IKaraali et al.l . 120071 ; iBilir et al. 



2008al ). The dependence of the local space densities of thick disc and halo 
on galactic latitude for high latitude fields is a result of a bias effect. High 
latitudes give the chance to the investigator to reach greater distances in the 
vertical direction, for a specific distance in the line of sight. That is, in high 
latitudes, the halo stars are dominant and they affect the galactic model as 
their contribution results in different local space densities, depending on the 
number of halo stars. 
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Fig. 17. Variation of the scaleheight of the thin (a) and thick (b) discs with R. The 
dots, the dashed and s o lid lin es represent estimated scaleheights, using adopting 
Lopez-Corredoira et al. ( 20021 ) 's flare model (Eq. 21) and calculated flare model 
using scaleheight of 20 star fields (Eq. 22). 

4-2 Scaleheight and Scalelength of the Thin Disc 



As cited in Section 3.3, the scaleheight and scalelength of the thin disc vary 
with galactic longitude, however, with different trend s. These variations are 
explained as a disc-flare. lLopez-Corredoira et al.l (120021 ) showed that the scale- 
height of the thin disc increases with galactocentric distance (R): 



h z = h Zj0 [l + 0.2l(R - Rq) + 0.056(i? - Rof 



(21) 



where Rq is the galactocentric distance of the Sun and h z>0 is the scaleheight of 
the thin disc at the solar distance. We plotted the scaleheights of the thin disc 
for 20 fields versus galactocentric distances in order to confirm this equation. 
Fig. 17 shows that the scaleheight of the thin disc increases with galactocentric 
distance and supplies the expected confirmation: 



#i = (349.68 ± 39.91)[1 + (0.053 ± 0.013)(i2 - R ). 



(22) 
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Fig. 18. Variation of the number of stars with go apparent magnitude intervals as a 
function of the galactic longitude. 



4-3 Scaleheight and Scalelength of the Thick Disc 



The variation of the scaleheight and scalelength of the thick disc with galactic 
longitude are due to: 1) the long bar, and 2) the flares. The bar induces a grav- 
itatio nal "wake" , traps and piles u p stars behind it (IHernquist &: Weinberg 



19921 ; iDebattista fc Sellwoodlll998l ). Thus, in response to a bar, one expects 
an excess of stars in the direction of the long bar (/ ~ 27° and ~ 207°) and a 
deficiency of stars in the opposite direction, which affects the numerical values 
of the scaleheight and scalelength of the thick disc. It is worth noting that the 
maximum surface density for stars with apparent magnitudes 15 < go < 18 
(where thick disc stars are dominated) corresponds to galactic longitude ~ 20° 
(Fig. 18). 



Also, the flare produces a change in the scaleheight of th e thick disc as the 



galactocentr ic dist ance varies. This effect i s well explained in lLopez-Corredoira et al 



( 120021 . 12004 ) and ICabrera-Lavers et all (120071 ) . who showed that the thick 
disc presents a flaring in the opposite sense of the thin disc. That is, the 
scaleheight of the thick disc decreases with inc r easing galactocentric distance 



( ILopez-Corredoira et all 120021 ; iMomany et all 120061 ) . We plotted the scale- 



heights versus the galactocentric distances of 20 fields and confirmed this 
suggestion. Actually, panel (b) in Fig. 17 shows that the scaleheight of the 
thick disc d ecrease s with galactocentric distance. We are quoting the work of 
Bilir et al. (l2008ah where the same confirmation was done. 
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4-4 Axial Ratio of the Halo 



The trend of the axial ratio (c/a) of the halo with galactic longitude is almost 
flat ((c/a)=0.56), i.e. its range for 75% of the fields is only 10°. However, there 
is a slight minimum and a small maximum corresponding to the fields in the 
second and third quadrants, respectively. The flatness of (c/a) can be exposed 
by the fact that the halo component of the Galaxy investigated in this work 
corresponds to the "inner halo". Since the maximum galactocentric distance 
of the fields is less t han the upper limit of the inner halo, i.e: R: 10-15 kpc 



( jCoroilo et al.l . 120071 ). this statement is confirmed. Also, the inner halo has 



a dislike structure which is confirmed by the mean axial ratio ((c/a) =0.56) 
estimated in our work. Additionally, the inner halo is the outermost part of the 
Galaxy which can be investigated by the use of intermediate galactic fields. 



The variation of the halo model parameters with galactic latitude and lon- 
gitude estimated in high latitude fields are usually explained by means of 
two c ompeting scenarios, as c i ted in the Sect i on 1: 1) the triax iality of the 
halo (INewberg fc Yannvl 120061 : IXu et al.l 120061: ljuric et al.ll2008h. an d 2) the 
remnants of some historical merger events (IWyse fc Gilmord 120051 ) . Fig. 18 
shows that the maximum and minimum surface densities for the apparent 
magnitudes 18 < go < 22, corresponding to halo component of the Galaxy, 
do not fit with the galactic longitudes / = 0° and I = 180°, respectively, con- 
firming these scenarios. Instead, the maximum surface density is at I = 20°. 
But, this value is close to the direction, which the long disc bar points out. 
This finding encouraged us to argue that the slight minimum and the small 
maximum cited above originate from the long bar effect. To confirm this ar- 
gument, one needs to remember that the maximum of (c/a) corresponds to 
the fields in the third quadrant, i.e. 210° < / < 220° where long bar lies, 
whereas its minimum is obtained for the star fields deficient in the second 
quadrant. Also, we should remind that we worked on intermediate latitude 
fields, where the halo is rather shallow. The m i nimum and maximum accre- 
tions, i.e. S agittarius streams (INewberg et al.l (120021 ). their Fig. 1) and/or 
Palo mar 5 (jOdenkirchen et all 120031 ) are centred at the galactic coordinates 
(0°,45°) and (20°, 45°), respectively. 



Three excesses of surface densities, (I = 20°), (I = 50°, 19 < g < 20), 
(I = 210°, 20 < go < 22), in Fig. 18 correspon d to the overden s ity re gions 
cited in the literature i.e. Herculus-Aquila Cloud (Belokurov et al. . 20071 ) . and 
Sagittarius Star Stream (INewberg fc Yannvl . 120061 ). Thus, the distribution of 
surface densities in Fig. 18 summarizes many properties of the galactic com- 
ponents. 
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4-5 Metallicity Gradient in the Vertical Direction 



The metallicity distribution could be obtained up to distances z = 10 kpc from 
the galactic plane, which covers the thin and thick discs and inner halo. There 
is an agreement between the metallicity gradients in the vertical direction 
evaluated in our work and the corresponding ones appearing in the literature. 
That is, for short distances, the metallicity gradient is d[M/H]/dz ~ —0.3 
dex kpc -1 , confirming the formation of the thin disc by dissipational collapse. 
The metallicity gradient is steeper in the transition regions of two galactic 
components, such as thin and thick discs, and finally the inner halo may be 
assumed as a component of zero metallicity gradient, which is subject to a 
different formation scenario, i.e. one expects a contamination from mergers or 
accretion of numerous fragments from objects such as dwarf galaxies. 
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